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I (54) Title: MULTI-PATH RESISTANT FREQUENCY -HOPPED SPREAD SPECTRUM MOBILE LOCATION SYSTEM 
(57) Abstract 

A multi-path resistant frequency-hopped spread spectrum 
mobile vehicle or personal location system is described which 
provides low cost manufacture and low power operation while 
still enabling the accurate location of the mobile unit over long 
distances and inmoderate to severe mutli-path conditions. The 
frequency-hopped spread spectrum mobile vehicle or personal 
location system consists of a cenral station, a plurality of base 
stations and a plurality of mobile transmitters which transmit 
using a frequency- bopped spread-spectrum differential bi phase 
shift keying communication signal. Frequency Shift Keying 
modulation may also be used. Each of the plurality of base 
stations include an array of receiving dipole antennas and employs 
special algorithms for retrieving very low power frequency- 
hopped spread spectrum signals in a noisy and multi-path 
environment. The base stations use computational algorithms for 
determining the phase difference between each of receiving dipole 
antennas to determine the direction of the trans miner relative to 
the location of the respective base station. The multiple direction 
of arrival angles of the received signal at each base station are 
corrected based on an n-dimensional ambiguity space to locate 
most probable angles of arrival. The ambiguity space plot is 
used to eliminate erroneous comparisons of dissimilar phases of 
the mcorning signal. The most probable values are then tracked 
over multiple frequency hops and a histogram analysis of the 
strongest surviving angles of arrival is performed. The two peaks J 
of the histogram are used as the two mots probable directions 1 10 

of arrival. Each base station then communicates the two relative 

directions of the transmitter to a central station where the location of the transmitter is determined by triangulation. The direction of amval 
angles from the multiple base stations are summed in a least mean square approach to find a single direction of amval which is then plotted 
on a gnomonic projection to correct for the curvature of the earth. 
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MULTI-PATH RESISTANT FREQUENCY -HOPPED SPREAD 
SPECTRUM MOBILE LOCATION SYSTEM 



FIELD OF THE INVENTION 
The presort invention relates generally to radio communication 
systems and in particular the present invention relates to vehicle or mobile 
locators using frequency hopped spread spectrum communication in a multi- 
10 path environment 

BACKGROUND OF THE INVENTION 
Location of radio transmitters is a technology that has been 
used for many years and in many applications. In general, radio direction 

15 finding is classified into many subclasses using a variety of direction finding 
tools such as antenna loops and antenna arrays. For radio location, the use of 
mobile transmitters and a plurality of fixed base station receivers is generally 
known in the art where the fixed receivers locate the source of the radio 
signal by triangulation of the radio signals. In the alternative, a plurality of 

20 fixed base stations may transmit radio signals and a mobile receiver may 
ctetermine its own position 

Triangulation is a technique where the fixed base station 
receivers each determine the approximate direction of the radio signal from 
the transmitter. The approximate direction lines are then drawn on a map 

25 where the lines will cross. Since the direction lines are approximate, and 
typically three fixed base station receivers are used, the intersecting area 
between the three lines is usually shaped like a triangle since the lines rarely 
meet at a single point The transmitter is then located within the triangle. 
Triangulation of radio signals is fraught with uncertainty. Propagation 

30 interference, signal drop out, insufficient transmitter power, ambiguity in 

bearing measurements received from behind the receiving antennas and multi- 
path reflections are common sources of errors and uncertainty in radio 
direction finding and location. 

Vehicle location systems are known in the art which use a 

35 variety of direction finding and location techniques. One example of a 
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vehicle location system is U.S. Patent No. 4^740,792 to Sagey et at. which 
describes a criohile transmitter and a plurality of satellite receivers. A central 
site detemrines the location of the mobile transmitters by calculating the 
difference in the arrival time of the spread spectrum transmitter signals at the 
5 satellites. This system may be used, for example, for the anti-theft location of 
stolen vehicles or the location individuals canning the transmitters. 

In the alternative, a plurality of fixed base stations may transmit 
radio signals and a mobile receiver may determine its own position as for 
example in Loran-C or Global Positioning Satellite Systems. One of these 

10 types of systems is exenplified in US: Patent No. 4,665,404 to Christy et al. 
which desoibes the use of spread spectrum radio transmitters at fixed 
locations and mobile receivers. The receivers locate their own position by 
accurately measuring the time at which the radio signals are received to 
determine the range from the transmitters. 

15 There is a need in the art for a ground-based vehicle or mobile 

transmitter location system which locates a large number of small, low-power 
transmitters and which operate without confusion or collision There is a 
further need for a ground-based mobile transmitter location system which is 
immune from multi-path disruption of the radio signals in an urban or 

20 irregular terrain environment There is a further need for a ground-based 
mobile transmitter location system which is capable of carrying messages or 
identification signals simultaneous with the location of the transmitter. 

For the reasons stated above, and for other reasons stated below 
which will become apparent to those skilled in the art upon reading and 

25 understanding the present specification, there is a need in the art for a 

low-power, spread spectrum vehicle or personal locator system which operates 
effectively in moderate to severe multi-path conditions without collisions of 
simultaneous messages from a plurality of transmitters. 



30 



SUMMARY OF THE INVENTION 
The above mentioned problems with vehicle or personal 
locators and other problems are addressed by the present invention and which 
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will be understood by reading and studying the following specification. The 
present invention is a ground-based vehicle, mobile or personal transmitter 
location system which locates any one of a large number of small, low-power, 
.frequency hopping spread spectrum transmitters and which operate without 
5 confusion or collision of the received signals. The present invention is 

- immune from multi-path disruption of the radio signals in an urban or 

irregular terrain environment is capable of carrying messages or identification 
signals simultaneous with the location of the transmitter. 

The present invention consists of a central station, a plurality of 

10 base stations and a plurality of mobile transmitters which transmit using a 
fiequOTcy-hc^ped spread-spectrum communication signal. The 
communication signal may be modulated using differential tri-phase shift 
keying, or frequency shift keying, or the like. Each of the plurality of base 
stations include an array of receiving dipole antennas and employs special 

15 algorithms for retrieving very low power frequency-hopped spread spectrum 
signals in a noisy and multi-path environment The base stations use 
computational algorithms for determining the phase difference between each 
of receiving dipole antennas to determine the direction of the transmitter 
relative to the location of the respective base station. The multiple direction 

20 of arrival angles of the received signal at each base station are corrected based 
on an n-dimensional ambiguity space to locate most probable angles of 
arrival. The ambiguity space plot is used to eliminate erroneous comparisons 
of dissimilar phases of the incoming signal. The most probable values are 
then tracked over multiple frequency hops and an histogram analysis of the 

25 strongest surviving angles of arrival is performed The two peaks of the 

histogram are used as the two most probable directions of arrival. Each base 
station thai communicates the two relative directions of the transmitter to a 
central station where the location of the transmitter is determined by 
triangularion. The direction of arrival angles from the multiple base stations 

30 are summed in a least mean square approach to find a single direction of 
arrival which is then plotted on a gnomonic projection to correct for the 
curvature of the earth. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
In the drawings, where like numoals refer to like components 
idirou^cni the several views, 

Figure 1 is an overview of the infrastructure of the multi-path 
5 resistant freqikiKy-hopped spread spectrum mobile location system. 
Figure 2 describes the synehrarizati 
the outgoing paging signals from the base stations. 

Figure 3 describes the format of the freqpjmcy-hoppbd spread 
spectrum signal transmitted by the remote mobile units. 
10 Figure 4 is a diagram of the receiving antennas of the base 

stations including a plurality of dipole antennas and a single reference omni 
r - antenna. 

Figure 5 is a three antenna array for an interferometric direction 
finding system 

15 Figure 6 is an ambiguity plane plot of the three; antenna array 

of Figure 5. 

Figure 7 is the ambiguity plane for the ratios 5:6 of 7J1 the 
three antenna array of Figure 5. 

Figure 8 is an example of multi-path reflections in an urban 

20 environment. 

Figure 9 is diagram showing triangulation of the location of a 
remote mobile unit between three base stations. 

Figure 10 is a diagram of a gnomonic plot of the measured 
bearings of the direction finding algorithm on a sphere. 

25 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
In the following detailed description of the preferred 
embodiment, reference is made to the accompanying drawings which form a 
part hereof, and in which is shown by way of illustration specific preferred 
30 embodiments in which the inventions may be practiced These embodiments 
are described in sufficient detail to enable those skilled in the art to practice 
the invention, and it is to be understood that other embodiments may be 
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utilized and that structural, logical and electrical changes may be made 
without departing from the spirit and scope of the present inventions. The 
following detailed description is, therefore, not to be taken in a limiting sense, 
and the scope of the present inventions is defined only by the appended 
5 claims. 

The present invention is directed to a location system for use in 
locating remote mobile units. Hie remote mobile units may be located in 
motor vehicles, located on the person of people, carried within containers or 
packages, or any number or variety of mobile carriers. The location system 

1G comprises a plurality of base stations (at least two) which transmit signals to 
the remote mobile omits and receive signals from the remote mobile units. 
The base stations may send messages and information to the remote mobile 
units, or the base stations may simply said alert or interrogation commands to 
activate the remote mobile units. Thus the remote mobile units may be 

15 activated locally (for example by the wearer or vehicle) or remotely by the 
base station. The remote mobile units may be used by a person as a reverse 
pager, as an emergency locator or as a communication device. The remote 
mobile units may also be used as part of a vehicle to locate the vehicle if 
stolen, in an accident, for vehicle tracking or as an integral communication 

20 device. 

The remote mobile units receive messages as a standard paging 
device over licensed airwaves using a standard p aging infrastructure. The 
remote mobile units, when activated, transmit low-power (less than one watt), 
frequency-hopped, spread-spectrum communication signals. The transmitted 

25 signals from the remote mobile units are received by the base stations which 
are thai used to locate the remote mobile units. The base stations are 
equipped with an array of irregularly-spaced dipole receiving antennas and a 
reference antenna (optional). The direction of the received signals from the 
remote mobile units is determined by interferometric direction finding using 

30 phase difference comparisons between the various dipoles of the array. 

Unique algorithms using histogram analysis and spatial filtering are performed 
on the received signals at the base stations to eliminate the ambiguity caused 
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by signal interference, signal drop out, ambiguity in beating measurements 
received from behind the receiving antennas and multi-path reflections. 



System Overview 

5 The infrastructure of this acknowledgement paging systems 

serves as the basis of the present invention and the reverse pagers of this 
acknowledgement pager system operate identically in the present system The 
present invention allows the accurate location of the reverse pagers even in 
noisy and multi-path environments such as those found in urban environments 

10 and irregular terrain environments. 

Figure 1 depicts the major components of the two-way paging 

■ system in copending US. Patent i^hcatim Number 08/158,441 entitled 
"LOW-POWER FREQUENCY -HOPPED SPREAD SPECTRUM 
ACKNOWLEXXSMENT PAGING SYSTEM' filed November 24, 1993. In 

15 the preferred embodiment of the present invention, all of the components of 
the existing paging system of Figure 1 are used The reverse paging terminal 
1 10 at the central site operates to provide synchronization and messaging 
information through the paging terminal 116 to the reverse pagers 100 (also 
known as remote mobile units 100) via direct links to the base stations BS,, 

20 BS 2 and BS3 through ground based radio links (not shown) or through a 
satellite uplink/downlink using a geostationary satellite (not drown). The 
base stations BS,, BS2 and BS3 include transmit and receive towers 1 13a, 
113b and 113c, respectively and base station terminals 201a, 201b and 201c, 
respectively. Terminals 201 a-201c are required for producing the accurate 
25 synchronization information needed to be transmitted to the remote mobile 
units and for local processing of the received messages for direction finding. 
This synchronization information is used to coordinate the response of 
messaging from the plurality remote mobile units 100 so as to minimize 
collisions within groups of remote mobile units and eliminate collisions 
30 between groups of remote mobile units. 



***** 
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Base to Remote Svn chrnniTatinn 

Standard paging messages seat from the base stations BS,, BS2 
and BS3 to the plurality of remote mobile units 100 are, in the prefeired 
embodiment, sent as digital data encoded in the POCSAG paging standard 
5 These messages may be? used to interrogate the remote mobile units to activate 
the remote mobile unit to allow the base stations to begin the location process. 
Typically the paging channel has a center frequency of 143. 160 MHz, with an 
NRZ FSK data rate of 5 12 bps or 1200 bps. Other bit rates such as 2400 
baud (bps) are also feasible. Figure 2 describes the POCSAG paging 

10 communications protocol as modified for use by the prefeired embodiments of 
the present invention. In the top line of Figure 2, a greatly compressed time 
line of digital data transmitted acceding to the POCSAG protocol is shown. 
Batches of messages are transmitted in groups as shown in the details in the 
subsequent lines below the top line of Figure 2. In the second line of Figure 

15 2, a 1.0625 second interval (for 512 baud) is shown in which 544 bits are 
transmitted as a single batch. The batch is preceded by a synchronization 
code word SC as shown in the third line of Figure 2. This synchronization 
code word is used to get the attention of remote mobile units 100 in the 
geographic locale serviced by the paging terminal. 

20 The synchronization code word within each batch is followed 

by eight frames of digital data. Each frame is divided into two portions, an 
address portion and a message portion. The address code word of the 
message of frame 2 of Figure 2 is shown in line 4 while the message code 
word of the second half of frame 2 is shown in line 5. The address code 

25 word is preceded by a digital zero followed by 18 address bits, two function 
bits and 10 check bits. The address code word is followed by an even parity 
bit The message code word portion of the frame is preceded by a digital one 
followed by 20 message bits which are followed by 10 check bits and a single 
even parity bit. Thus each frame is comprised of 64 bits divided into two 32 

30 bit sections. 

Synchronization of the base station terminal 200 and the remote 
mobile units 100 is necessary to ensure the units 100 are transmitting at the 
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same time that the bas e ^ Synchronization is also 

necessary to coordinate the division of the - large number of remote mobile 
units into groups so that members of one group use different frequency 
helping p Syndnabi^ the 

5 %edK^ffi6fcfe^ inserting a special frame info the 

AG data t^ci i^ 

Synchronization brtweai the reverse paging terminal 110 and 
the remote mobile units 100 is accomplished on tw^ Synchronization 
of the remote mobile units within groups of remote mobile units determine 

10 vvhere along the pseudo random noise code the frequency hops are to be 
followed. For example, within a singje group of remote mobile units, all of 
the remote mobile units within that group will be synchronized to begin 
transmitting at the same location in the pseudo random noise code list for any 
acknowledgement which may be required Synchronization information is 

15 sent from the reverse paging terminal periodically to the addresses of each of 
the remote mobile units within each group to remind the remote mobile units 
100 where along the pseudo random noise code they should be tracking. This 
also enables the dynamic changing of a remote mobile units group 
membership such that if one group is experiencing a large number of 

20 collisions due to simultaneous transmissions, the reverse paging terminal 110 
may re-allocate some of the remote mobile units within that group to new 
groups to minimize collisions. 

Another form of synchronization is required to synchronize the 
remote mobile units to the exact times for transmitting frequmties from 

25 within any of the hops. This fine synchronization information, described 
more fully below, is transmitted as part of the POCSAG codes. 

Referring once again to Figure 2, eight frames of information 
are transmitted in each burst using the POCSAG format. Remote mobile units 
100 may be assigned to a specific frame within the transmission so that the 

30 remote mobile units, once recognizing the synchronization code word, can 
scan a specific frame for that remote mobile unit's address. Once the address 
is found, the remote mobile unit can determine any group changes that may 
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be required to re-allocate that remote mobile unit to a different group. In 
addition, the POCS AG format is used to transmit a fine time synchronization 
code. The fine synchronization code is a transmission of a time pulse at an 
exact time syntonized to a GPS (Global Positioning System) clock to 
5 syndironizB all the remote mobile units 100 for time of transmission For 
example, periodically during the day the reverse paging temrinal will send a 
syndarcMiization code within the POCS AG code wond winch is sent at a very 
precise time: In order to ensure that a precise time pulse is sent, the reverse 
paging terminal 110 receives accurate time information using a GPS antenna 
JO to receive accurate time of day information. The time used to send the 

synchronization pulse is when the day clock reaches exactly some multiple of 
0.9 seconds in the prefered embodiment. In this synchronization information, 
20 bits of information are transmitted to give the accurate time of day 
information. . - ■ 

15 In each of the remote mobile units 100, the microprocessor 

compares this accurate time pulse which will indicate the exact time of day 
and compare it to its own day clock. The clock within each microprocessor is 
accurate down to a few milliseconds, but the time at which the 
synchronization pulse occurs should have a resolution much finer than that 

20 such as down to 0.1 milliseconds for time of day. In this fashion, each of the 
microprocessors in each of the remote reverse paging devices can periodically 
realign its day clock to know within a millisecond the exact time. Each 
microprocessor does not actually realign its clock but changes a clock offset 
within memory so that it understands how far off its own internal clock is and 

25 can make the adjustment when using that clock to determine when to start 
transmitting information using the eight frequency hopped spread spectrum 
signal. 

The synchronization pulse is only transmitted every few 
minutes. However, the resolution of the start of the message indicating the 
30 synchronization pulse is very accurate, it being transmitted at 0.090000 

seconds GPS time after a fixed time of day, such as 12:00 GMT. This GPS 
time is accurate to at least within 100 nanoseconds. 
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An overview of the transmission fonnat of the remote mobile 
unit is shown in Figure 3. The actual transmission of information from the 
remote mobile units 100 is done ming Differential Bi-Fhase Shift Keying 
(DBPSK) modulation on a frequency hc^^ 
5 The transmission of information from the remote mobile units 100 on the 
frequency hopped carrier may also be done using Frequency Shift Keying 
(FSK) modulation. Typically a single transmission ^ consists of 53 hops or 53 
changed frequencies selected from a list of narrow band frequencies. The 
frequency selection is based on a pseudo-random noise code list pointing to 

10 the frequency selection list The synchronization infoimatian tells the remote 
mobile unit 100 where along the pseudo random noise code it should be 
synchronized for transmission of its message and the fine synchronization 
information tells exactly the time of day so that the remote mobile unit 100 
knows exactly when to start transmitting the specific frequency so that die 

15 reverse paging temrinal 1 10 is looking for that frequency at die same time. 

In operation, 200 frequencies are used by the remote mobile 
unit 100 and the base stations and internally stored in a list numbered Fl 
through F200. For a specific message, 53 frequencies will be used to transmit 
the entire message. These 53 frequencies are selected based on a 1 ,000 

20 member pseudo-random noise code. 

The use of the accurate syndhronization signal periodically 
broadcast via the outbound paging signal enables the remote mobile units to 
use lower accuracy components thus reducing the manufacturing cost of 
remote mobile units. For example, high accuracy crystals to track the time of 

25 day within the microprocessor are available with an accuracy of three parts 
per million. Thus, a time drift of approximately three micro seconds per 
second or 180 microseconds in a minute is the known drift. There are also 
time inaccuracies which are introduced due to the time of transmission 
(variable path length) from die source from the reverse paging terminal when 

30 the synchronization information is sent. By employing crystals which are 

cheaper and have an accuracy of the order 50 parts per million, the amount of 
time-of-day drift normally wouldn't be tolerable. However, by using the 
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synchronization information transmitted on a regular basis from the reverse 
paging terminal, the microprocessor can continually correct its own internal 
day clock so that accurate time of day measurements are always maintained 
The microprocessor estimates the momentary inaccuracy of the crystal by 
5 t racking the drift across several synchronization transmissions and 

dynamically adjusts for the frequency drift of the crystal and the offset using 
internal offset registers for accurate time of day information. 

Counters are employed within each microprocessor of the 
reverse paging units to compensate for the oflset of the frequency based an 

10 the synchronization time information There are generally two major factors 
which affect the drift in a oystal: temperature and acceleration. Most of the 
drift is due to temperature, and the remaining drift components are negligible. 
The frequency drift in a oystal due to temperature is very slow, on the order 
of 50 Hz over 10 seconds. During a single day die temperature can change 

15 by 20 or 30 degrees fahrenheit, requiring a time update from the GPS clock 
approximately every five minutes. 

Remote Mobile Unit Transmission Format 
The signal sent from the remote mobile unit 100 to the base 

20 stations is a spread-spectrum, frequency-hopped transmission using differential 
bi-phase shift keying (DBPSK) modulation on the frequency-hopped carrier to 
transmit digital information- The transmission of information from the remote 
mobile units 100 on the frequency hopped carrier may also be done using 
Frequency Shift Keying (FSK) modulatioa The frequency hops are relatively 

25 slow, the frequencies transmitted are very narrow and the transmission power 
is extremely small. The maximum peak output power of transmission from 
remote mobile unit 100 is limited to less than one watt to allow use of the 
902-928 MHz band in the United States without the need for licensing the 
paging transmitters as allowed by FCC regulations defined in 47 C.F.R 

30 §15.247, which is hereby incorporated by reference Those skilled in the art 
will readily recognize that other frequency bands and transmissions power 
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levels may be employed depending upon FCC licensing requirements or other 
frequency licensing reqtdrements of other nationalities. 

Ibe use of an accurate crystal to control each frequency of 
trahsmk^ 100. For example, 

5 high accuracy crystals to transmit the narrow bandwidth frequencies used for 
the frequency hopped transmissions are available with an accuracy of three 
parts per milli on At 900 MHz, a 3 ppm drift would place a single fiequency 
somewhere within a 2.7 KHz band. To tolerate fiequency drift due to aging 
and teb^)erature, each individual frequency of the frequency hopped signal is 

10 allocated to a 7.5 KHz band or channel, even though the actual frequency is 
on the order of 200 Hz wide skirt within this 7.5 KHz allocated bandwidth. 
Those skilled in the art will readfly recognize that by us^ 
components, the frequency channels (individual frequency of the frequency 
hopped signals) of 7.5 KHz allocated bandwidth may be wider or more 

15 narrow depending upon the overall allocated bandwidth for the system. For 
example, 1 KHz or less bands may alternatively be allocated per channel 
Tests on this invention have shown that by processing the 
received signals at the base stations entirely in the digital domain using the 
combination of unique Fast Fourier Transform algorithms of the present 

20 invention to locate and retrieve the frequency hops and by using a 
combination of unique confidence algorithms with a plurality of error 
correction codes, the receiving base station is able to pull the response 
information from a very low power signal from a distance of up to 45 
kilometers (28 miles) in a flat terrain- In a rather noisy urban environment, a 

25 range of 24 kilometers (15 miles) is the norm The information within the 
signals is accurately decoded even though direction of the incoming signals 
may be severely distorted due to multi-path reflections and noise 

As shown in Table 1, the remote mobile unit message format 
consists of a preamble and the message body spanning a total of 53 frequency 

30 hops. Those skilled in the art will readily recognize that longer messages may 
be transmitted using the preferred embodiment of the present invention, and 
the messages format described here is illustrative but not limiting. Since FCC 
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regulations defined in 47 CER §16.247 require a minimum of 50 frequency 
hops, the 53 hop message format described here is considered a short message 
hop. Much longer message h^ 

implemented but not described here. Of course, those skilled in the art will 
5 readily recognize that shorter messages than those described below are equally 
possible for the preferred embodiments of the present invention The message 
length and number of transmission hops are a matter of design choice. 

The message preamble consists of a predefined code of ones 
and zeros to get the attention of the base unit receiver to begin its FFT (Fast 
10 Fourier Transform) routines to begin pulling the message out of the noise. 
The preamble consists of 165 bits transmitted across 5 hops, that is, 
(\ / transmitted using DBPSK (Differential Bi-Phase Shift Keying) cm* Frequency 

Shift Keying (FSK) on five different frequencies selected from the frequency 
list with the specific frequencies selected based on the FN (Pseudo-random 
15 Noise) Code list stored within the remote mobile unit The sequence location 
within the PN code that the remote mobile unit will begin to follow is based 
on the synchronized time of day. Within a single hop (a single carrier 
frequency), the carrier phase is modulated 33 times to encode the predefined 
one-zero pattern of the preamble. 
20 The message body follows the preamble and consists of three 

groups of data Each group consists of 30 actual data bits so that the entire 
message is, in the preferred embodiment of the present invention, 90 total data 
bits (although other bit length messages may be chosen). The actual data 
encoded within these 90 bits is described above and may be in any convenient 
25 coded format. Those skilled in the art will readily recognize that a wide 

variety of message formats and encoding of the data bits may be used without 
departing from the spirit and scope of the present invention. The encoding 
described here, however, has been proven effective in retrieving the data bits 
buried in background noise with a high degree of accuracy and a low actual 
30 error rate. 
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TABLE 1 : Renxrte Mohile Unit Message Fonmt 

Preamble is 165 bits (33 bits x 5 hops) 
Message is 512 transmitted bits 
(Message is 90 bits actual data) 

/7 = One Frequency Hop 
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Outer Message Coding 
Each of the three groups of message data (30 bits each) are 
15 BCH encoded using a standard 30,63 BCH code and with a single parity bit 
added to form a 64-bit word This encoding decreases the eiror rate from 10 2 
to 10" 5 . This encoding, documented and understood by those skilled in the art, 
can correct 15? to 6 errors or detect up to 13 errors. Detection of corruption of 
a data word that cannot be reconstructed will cause die base to request a 
20 second transmission of die acknowledgement message. 

Inner Coding and Int erleaving 
The inner coding of the message will protect the integrity of the 
message with an error rate as high as 25% Each block of 64 bits of data 

25 (corresponding to a groups of 30 bits and earlier encoded by a standard 30,63 
BCH code) is split into two sub-blocks of 32 bits (sub-blocks A and B of 
Table 2), and a reference bit is added to each sub-block to assist the 
differential encoding to provide a reference bit to the DBPSK or FSK 
decoder. The 33 bit sub-blocks are transmitted over one frequency hop each 

30 and are replicated 8 times so that the 64-bit block traverses 16 frequency 
hops. In transmission, the 33 bit sub-blocks are interleaved to farther reduce 
loss of data, as shown in Table 3, where sub-blocks A and B of Table 2 
correspond to the first group of 30 bits, sub-blocks C and D, correspond to the 
second group of 30 bits, etc. The total message is 53 hops where each hop is 

35 ] 80 msec in length making the duration of a single message 9.54 seconds. 
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TABLE 2: Interleaving Format for Sub-Hock 

A = 1 reference bit and 32 data bits = 33 bits 
B = 1 reference bit and 32 data bits = 33 bits 

[] = One Frequency Hop 
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TABLE 3: Inner Coding and Interieaving of Sub-Hocks 

A — first 33 bits of 1st block 
B = second 33 hits of 1st block 
C = first 33 bits of 2nd block 
D = second 33 bits of 2nd block 
E = first 33 bits of 3rd block 
F = second 33 bits of 3rd block 

[] = One Frequency Hop 
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Those skilled in the art will readily recognize that a wide 
variety of data interleaving may be utilized to effect better error tolerance and 
may be substituted for the interleaving described here. Such alternate 
substitute interleaving means are QRC (Cross Interleaved Reed Solomon 
Code) used in CD (Compact Disc) recording media operating either at the 
block level or at the bit level 



Single Hop Format 
The acknowledgment signals are transmitted by the remote 
35 mobile units 100 in a 1.5 MHz band selected from within the 902-928 MHz 
spectrum The 1.5 Mhz band is divided into 7.5 KHz channels to provide 200 
channels available in which the frequency hops can occur. Thus, each 
frequency hop is a channel 7.5 KHz wide in which a carrier frequency is 
transmitted. For example, channel one will have a frequency Fl at 
40 902.000000 MHz +/- 3.25 KHz, channel two will have its center carrier 
frequency at 902.007500 MHz +/- 3 KHz, etc. 
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Each transmit fr&qaeoty 
the approximate mid-point of the assigned channel band; however, due to 
inaccuracies in the remote mobile unit circuits and reference crystals, the 
actual transmit frequencies will Vary between units. If high quality crystals 
5 are used to accurately produce the required frequencies, very little drift off the 
cento- frequency will result In the preferred embodiment of the present 
invention, low cost crystals are purposely employed to keep the per-unit 
manufacturing costs down: This will allow few: a lower-cost product sold to 
the user which will increase market penetration. Thus, reference crystals are 

1 0 preferred which have a frequency accuracy of 3 ppm such that at 900 MHz, 
the statistical drift would be approximately 2700 Hz. The crystals center 
frequency within its nominal accuracy also drifts due to aging and temperature 
variations, but this drift is slow compared to the transmission times so the 
drift during a single transmission due to these latter variants is unimportant 

15 A single frequency hop is shown in Table 4. The 15 

millisecond guard time preceding each hop is primarily a settling time for the 
oscillator circuits of the remote mobile units to allow the internal oscillator 
circuit to lode onto the new frequency between hops. Each hop is transmitted 
at a single frequency in which the phase of die carrier is either at 0 degrees 

20 phase or 180 degrees phase in reference to the phase of the reference bit 
immediately following the quiet or guard time. Thus the first bit is a phase 
reference bit followed by 32 data bits exhibiting either zero phase shift or 180 
degree phase shift to encode the data bits as DBPSK (Differential Bi-Phase 
Shift Keying). In an alternative implementation, each frequency hop may be 

25 modulated using Frequency Shift Keying (FSK) in which two frequencies are 
used to transmit data bits. One hop frequency may indicate a logical one 
while a second hop frequency may indicate a logical zero. The frequency 
shift is minor and the frequency differential is contained within a single hop 
channel. 

30 Each bit of DBPSK or FSK is a transmission of approximately 

5 milliseconds of the hop carrier frequency either in phase with the reference 
bit transmission or 180 degrees out of phase. Each actual bit is approximately 
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4:7 milliseconds of carrier at the hop frequency preceded and followed by 
approximately 0.15 milliseconds of quiet guard band to reduce discontinuities 
between phase changes which contribute to step-function noise in the 
transmission. 

5 TABLE 4: Single Frequency Hop Format 

Guard time (quiet) = 15 ms 

Single Bit = 5 ms of carrier DBPSK / FSK 

33 Bits plus guard time = 180 ms 

to 
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Frequency Hopping Sequence 
All of the remote mobile units in the market serviced by the 
reverse paging terminal for message or location finding use the same pseudo 

20 random noise code to determine the frequency hops. The pseudo random 
noise code is a digital code which is repeating after approximately 1,000 
unique codes. In the preferred embodiment of the present invention, the 
pseudo random noise code is stored in memory of each of the remote mobile 
units. Those skilled in the art will readily recognize, however, that a linear 

25 feedback shift register could be used to generate the pseudo random noise 
code on a real-time basis instead of using a look-up table which is presently 
in the preferred embodiment. 

The PN (pseudo-random noise) code list is stored in memory 
and maps to a frequency list In the preferred embodiment of the present 

30 invention, the PN code list has 1,000 entries which repeat as a sequence. The 
control means of the reverse paging units continuously maintain a count of the 
proper location within this list for the exact time of day. As described below, 
the time of day for all remote mobile units in the market served by the base 
terminal are periodically synchronized to ensure acknowledgment messages 

35 are synchronized to transmit the hop frequency at the proper time and to 

synchronize the location within the PN code list that each remote mobile unit 
will use to transmit. 
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The 1,000 member PN code list maps to a 200 ihember 
frequency list In order to aUow a large number 
simultaneously operate in the same geographic market, the remote mobile 
units are divided into groups and the groups are assigned different sequence 
5 segment locations iti the same 1,000 member FN list Thus a remote mobile 
unit from group one will begin transmitting a hop at a frequency determined 
from a first location with the PN code, while a remote mobile unit from group 
two may begin transmitting a hop at a frequency detenu^ 
location in the PN code. The remote mobile units from group one and group 

10 two will complete their respective acknowledgement messages in 53 hops. 
Preferably, the sequence of the PN code used to determine the frequencies of 
the 53 hops for the remote mobile unit of the first group will not overlap the 
sequence of the PN code used to determine the frequencies of the 53 hops for 
the remote mobile unit of the second group. More preferable, the frequencies 

15 chosen based on the non-overiapping segments of the PN code list are 

orthogonal such that the same frequency is never used by two remote mobile 
units belonging to different groups. 

In the preferred implementation, the 1,000 member PN code list 
is divided into 160 hopping sequences. The remote paging units are divided 

20 into 40 groups with the members of each group synchronized to track the 

same location in the PN code list. The microcontroller of each remote mobile 
unit, regardless of its group membership, continuously runs through the 
repeating PN code sequence to stay in synchronization with the base unit and 
all other remote mobile units. Each group of remote mobile units is further 

25 divided into four subgroups such that the remote mobile units within each 
subgroup are assigned one sequence within the PN code list Although the 53 
hop sequence needed for each acknowledgement transmission may overlap the 
53 hop sequence used by a remote mobile unit in another subgroup., the 
transmission sequences of a remote mobile unit of one group is chosen to not 

30 overlap the 53 hop sequence used by a remote mobile unit in another group. 
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Base Station Design 
As described in the parent patent application, the analysis and 
decoding of the signals received by the base stations from the remote mobile 
units is done almost entirely in the digital domain. The earner frequencies of 
5 the frequency hops are down-converted to a lower frequency in each base 
station and are then digitally sampled The digital samples are? then processed 
to locate the phase information of interest for direction finding and message 
decoding. Each base station is constructed with a plurality of digital signal 
processor pipelines which enable simultaneous message decocting and 
10 direction finding of a plurality of simultaneously transmitting remote mobile 
units. 

As shown in Figure 4, the receiving antennas of the base 
stations are constructed using a plurality of dipole antennas 401a, 401b, 
though 401n (generally referred to as 401) arranged in an array, as shown in 

15 Figure 4 and a single reference omni antenna 404. The omni antenna 404 is 
not necessary for the preferred embodiment of the present invention since any 
one of the dipole antennas of Figure 4 may act as a reference antenna for any 
of the other antennas. The omni reference antenna 404 is used as a reference 
antenna because the omni antenna 404 exists as part of the reverse paging 

20 infrastructure upon which the present invention relies and is used for receiving 
the messages from the remote mobile units. 

In the preferred embodiment of the present invention, the 
direction finding (DF) antennas of each base station uses twenty-four sector 
direction finding dipole antennas 401 arranged in four linear arrays as shown 

25 diagrammatically in Figure 4. Each of the four linear arrays consists of six 
dipole antennas. The phase between the antennas is measured indirectly - 
each instant of time one of the 24 antennas is selected and its phase is 
compared to the omni antenna. 

The switching of the antennas is done every one bit (5 

30 milliseconds) of every frequency hop and is synchronized to the incoming 

data The switching occurs at the 0.7 millisecond gap between the transmission 
of the bits. The antennas 401 are continuously scanned since there could be 
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many remote mobile units transmitting simultaneously from different 
directions. Differentiating between the various simultaneously transmitting 
remote mobile units is capable. If the simultaneously transmitting ranote 
mobile units belong to different groups, they will be using orthogonal 
5 frequencies because they are using different sequences in the PN code. If the 
simultaneously transmitting ranote mobile units belong to the same group, 
they are nonetheless distinguishable due to the partic^ar xmque frequoocy 
offsets caused by the frequency drift off center frequency of the transmitter 
reference crystal. As described in the above referenced parent application, the 

10 simultaneously received signals may be transmitted within the same 7.5 KHz 
channel during a hopping sequence of remote mobile units from the same 
gxstsp, but the base stations are capable of distinguishing and tracking the 
hops of the simultaneously transmitting ranote mobile units. Ibis is due to 
the resolution of the Fast Fourier Transforms (200 Hz bins) used in the 

15 pipelined receivers, the narrow frequency skirt of the transmitters in the 
remote mobile units and the drift tracking performed in the base stations 
which can distinguish the unique transmission signatures of the simultaneously 
transmitting remote mobile units. 

The direction finding process, described more fully below, 

20 compare the exact phase amplitudes of each of the dipole antennas of the four 
arrays to the phase amplitude at the reference antenna In the antenna array 
shown in Figure 4, a circle surrounding the array is divided into eight sectors. 
The actual direction computation is done in the sectors according to Table 5. 

25 

TABLE 5: Sector Allocation on the DF Antenna Array 

30° to 60° - top and left arrays 

60° to 120° - left array 
30 120° to 150° - left and bottom array 

1 50° to 210° - bottom array 

210° to 240 b - bottom and right array 

240° to 300° - right array 

300° to 330° - right and top array 
35 330° to 30° - top array 
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Those skilled in the art will readfly recognize that tte 
direction finding array does not have to be installed in a single location Few- 
instance, each linear array of six dipole antennas could be installed on a 
different face of roof of a high rise buflding. In some cases a remote mobile 
5 unit cannot transmit from specific sector. For instance one of the sectors 
feces the sea In that case the direction finding array could be reduced, and 
only a single direction finding antenna will be instaUed at that sector. Usually 
the interfearometer does not operate in case that a remote mobile unit is very 
close to the base station (less than 1 kilometer). In such a condition, the 

10 other base stations will still perform the direction finding and the central site 
will still be able to perform the triangulation 

For the reasons exemplified more fully below, the spacing of 
the six dipole antennas in each of the four sides of the array is critical to 
maintaining a low degree of ambiguity in direction finding and multi-path 

15 rejection The base ratios of the base lines was found by simulations with 
several relative prime numbers and is selected to be 6 : 10 : 4 : 3 : 5 to 
produce a phase gain PG = 567L The ratio was computed through extensive 
simulations. Although the spacing between the various antenna is critical to 
accurate determination of the incident angle of the carrier, the order of the 

20 spacing of the antennas in the array is not critical. The order chosen here was 
based primarily upon a physical implementation consideration that the array 
frame be mechanically balanced 

Merferometric Direction Finding 
25 The direction finding technique used in the preferred 

embodiment of the present invention is performed by measuring the electrical 
phase difference between the antenna of the array caused by the angle of the 
incoming wavefront of the spread spectrum frequency hopped sig n al . In ideal 
conditions where there is no multi-path disruptions of the wavefront, the 
30 incoming wave from the remote mobile unit is coming in angle of 0 from the 
perpendicular line to the base line on which the antennas are mounted 
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Figure 5 is an example of how the interferometric direction 
finding technique of the preferred embodiment of the present invention is 
practiced Figure 5 shows, by way of example and not by limitation, only 
three antennas of the array. The three antennas are used to determine phase 
5 differences of the incoming signals between pa^ 

described above, the phase difference may also be measured between each 
antenna and a reference antenna such as the omni antenna 

The interimraeter direction finding (DF) technique calculates 
the direction of arrival of the wavefront by measuring the phase of the 

10 incoming wave front in several places in space. This interferometric method 
is a planar interferametric method which assumes that the antenna array and 
the transmitting remote mobile unit are all in the same relative plane. 
Differences in height of the remote mobile unit and the antenna array are 
ignore and cause little error in the direction finding of the present invention 

15 As seen from Figure 5, the three antenna 501, 502 and 503 are 

irregularly spaced such that distance between antennas 501 and 502 is 
great than distance D 12 between antennas 502 and 503. The distance between 
the antennas is selected to be a multiple of a prime number sequence to 
ensure an irregular spacing of all antennas in the array which increases the 

20 probability that an electrical phase difference will always be measured 

between some of the antenna members of the array. The distance between the 
antennas is important to minimize ambiguity errors in measuring the phase of 
the incoming wavefront If the antennas are spaced too widely apart, there is 
increased ambiguity in determining whether the phase difference of the same 

25 wavelength is being compared or whether the phase of different wavelengths 
are being mistakenly compared Since the frequency hopped spread spectrum 
signal of the present invention is selected to operate in the 900 MHz band, the 
wavelength A, of the incoming signal is approximately 30 centimeters. 

Again, referring to Figure 5, the phase difference between the 

30 two antennas 502 and 503 is know as D 12 and the formula for measuring the 
angle of the incoming signal is : 
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(1) (<p 12 +2lliq 2 ) /2*=D 12 SIN(e)/\ 

:•• where 

0 is the angle of the incoming wavefront measured from the 
5 perpendicular of the antenna array baseline, 

<ft 2 is the electrical phase measured between antennas 502 and 
503 

K 12 is an ambiguity factor since the system is capable of 

measuring (p between ± % since the antennas could be 
10 measuring the phase difference of more than one cycle 

D n is the distance between the antennas 502 and 503, and 
X is the wavelength of the incoming signal. 

According to this formula, as long as E\ 2 is smaller than A/2, the ambiguity 
15 factor K )2 equals 0. Thus the largest antenna spacing in the array is selected 

to be smaller that 15 centimeters. 

In order to estimate the inherent instrumental error of the 

direction finding interferometer, the above equation can be written in a 

different form: 

20 

( 2 ) SIN(Q) =((p 12 +27iiq 2 ) / (2uD x2 /A) 

Assuming that all the errors in estimating 0 is due to inaccuracies in 
estimating (p, differentiating the above equation produces: 

25 

(3) aid) =o{(p 12 ) / (2nD x2 COS(0) /X) 

where or is the standard deviation 
From this equation, inherent limitations in the system are determined. The 
30 electrical error is divided by the factor 271 E^ 2 / X which is termed the phase 
gain (PG) which reduces the error by the ratio between the distance and the 
wave length. Also, the error rapidly grows to infinity as 0 approaches 90°. 
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Thus the coverage of the present linear interferometer is limited to ± 60° as 
shown by the division of die array into sectors described in above in 
conjunction with Table 5. In order to reduce the error, the phase gain must 
be made as large as possible and each antenna face of the antenna array of the 
5 present linear intefeon^er array covers only 1/4 of the space. The antenna 
spacing in the array must have a distance smaller than 7J1 in order to resolve 
ambiguities. 

Ambiguities in Interferometer Direction Finding 

10 There are two types of ambiguities in interferometer DF: front 

to back ambiguities (the linear interferometer only covers effectively 1/3 of its 
surrounding area) and multiple possible solutions of the interferometer 
equations, in case where the distance between two antennas is greater than 
Wl. This latter ambiguity is caused by measuring the phase difference 

15 between different cycles of the incoming wavefront. 

The front to bade ambiguity is solved in the present invention 
by using multiple arrays, as shown in Figure 4. In order to cover 360°, the 
present invention uses four sub-arrays of six dipole antennas each, for 
coveting 90° The antennas in each array have a vertical beam width of 120° 

20 The base station scans all the antennas of the array through the 1:24 antenna 
selector 402 and selects the proper sub-array by comparing the amplitudes 
received from antennas in each array. If two arrays are almost equal in 
amplitude (within 3 dB), then the interferometer process is done on both 
arrays and the direction of arrival (DOA) of the signal is computed either 

25 from both or from the proper array which is selected according to the DOA 
closer to its perpendicular. 

The multiple cycles ambiguity is solved by using more than 
two antennas in each linear array By way of ©cample, but not by limitation, 
a three antenna array is as shown in Figure 5. The distances between the 

30 antennas that are elected as D^2*A/2 and D23=3*?y/2. Using these ratios, the 
equation labeled (2) above becomes: 

(4) SIN(Q) =(f> 12 /2n+K 12 
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(5) SIN(O) -9 23 / 37t+2iq 3 /3 

Figure 6 is an ambiguity plane plot of the three antenna array 
of Figure 5. Figure 6 shows the phase difference (p 12 plotted against the 
5 phase difference (p^ as an interferometer ambiguity plane plot Figure 6 

shows a family of lines* 601, 602, 603, 604, and 605, parameterized by K 12 & 
K23 with the same slope created as 0 moves from 0° to ±90°. The following 
formula are used to create the plane plot of Figure 6: 

10 (q>23/3*) +2^/3= (q> x2 /27i) +K 12 

The ambiguity resolution algorithm is as follows. The phase 
difference between the antennas is measured as (p 12 and (p^. On the 
interferometer ambiguity plane of Figure 6, find the line represented by K 12 

15 and K23 (ambiguity factors) that has the minimal distance, which is the actual 
probability distance, to 9^ and (p^. After solving die ambiguity, compute 
the exact DOA using the full span of the array i.e. 2.5A, which gives a phase 
gain(PG) of 5n. 

One could select ratios of N:N + l of A/2, where N is the 

20 number of antennas in the array This ratio of N:N + 1 of A/2 may be 

selected to be as large as possible in order to improve the phase gain. But as 
N gets larger, the distance between the ambiguity lines gets smaller and the 
algorithm becomes vulnerable to an error in the proper selection of K 12 and 
K zv The result will be a gross error. 

25 If the measurements distribution of (p is gaussian, then in order 

to compute the gross error probability, the "tail" of the gaussian distribution 
from the distance between the ambiguity lines to "infinity" must be computed 
since the phase is folding after ± n . This computation described above is 
only an approximation and is almost true for small values of a ((f)) which is 

30 true at high signal to noise ratio (SNR) or high signal to interference ratio. 
For white gaussian noise, a ((p) could be approximated by 1/VsNR where 
SNR is the signal to noise ratio. 
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Figure 7 is the ambiguity plane for the ratios 5:6 of A/2. The 
distance between the ambiguity lines became less than half compared to the 
ratio of 2:3 which practically increased the gross error to be impractical. 
Thus, the solution is to add mare antennas to the array. Adding more 

5 antennas at different ratios spreads out the ambiguity lines in hyper space 
(multi-^limeasional space) with dimension of N-I (N - number of antennas in 
the array) and increases the distance between them. In the preferred 
embodiment of the present invention, it was concluded that in an urban 
environment six antennas are required in a linear array resulting in a five 

10 dimensional plot of Figure 6 

drawing). The following formula are used to create the plane plot of Figure 7 
which is an example of moving the antenna spacing closer together: 

SIN(d) = a/Di 2 ) * (<p 12 +2*Je L2 )/ (2ii) = (<P 12 /5tO +2K 12 /5 

15 

SIN{d) = (A/P 23 ) * (<p 33 + 27Ci^ 3 ) / (2tt) = (<p 23 /6*) +i^ 3 /3 
(<p 23 /6tt) +iC 23 /3 = ((p 12 /57l)+2iT 12 /5 

Multi-Path Error Sources in Interferometric Direction Finding 
20 The preferred embodiment of the present invention provides 

accurate direction finding and location in urban areas. The major source of 
error is multipath distortion in which the incoming wave front is impaired by 
reflections coming from different directions. Figure 8 is an example of multi- 
path reflections in an urban environment. The unique feature of frequency 
25 hopping enables the interferometer system of the present invention to 

distinguish between the principal reflection and other reflections. Since it is 
not guaranteed that a line of sight exist between the remote mobile unit and 
the base station due to multi-path reflections, die most consistent reflection 
through the hopping message is the one closest to the real bearing. Empirical 
30 testing and observations proved that the bearing error is in practice the size of 
the "urban block" from which the remote mobile unit transmits. In practical 
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tests in an urban environment, the urban block size was about 75 meters in 
average, and up to 200 meters peak error. 

Multi-Path Resistant Direction Finding Algorithm 
5 The full direction fin ding algorithm relates to measurements of 

a linear antennas array having six dipole antennas per side of a four-sided 
array with a single omni reference antenna (see Figure 4). The calculations 
are performed on a six antenna subarray after the sub-array is selected to be 
closest to the remote mobile unit as described above in conjunction with 
10 Table 5. The full interferometer formula for a six antennas sub-array is: 
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where 

q>i is the electrical phase measured between the antennas, 

Dj is the distances between the antennas 

X is the wavelength 
20 9 is the direction of arrival (DOA) 

Kj is the ambiguity factor, integer value such that, 
-7C<27iD i * SIN(0)/X+27cK i <7t 

Ej is the additive noise (assumed to be gaussian). 

To solve the equation is (6) above, linear array calculations are 
25 performed to solve the following vector values: 
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The results of these calculations are as follows: 
20 the vector elements (p are limited by: I -0.5< q)j < 0.5 | ; 

g) is the measurements vector and d is known; 
-0.5<X<0.5; X is unknown; 

k is unknown; its elements are integer, limited by the elements 

ofd 

25 E is die noise vector that is unknown. 

In the calculations, it is assumed that the noise is additive such as: 
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(8) 



Eiz d }= 0; i e(l # 2,3,4,5) 
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Where : 

is the electrical phase measured at antenna i. 

<r =<!>»- 

The solution to this equation is to search among all K e N, the 
20 specific K that minimizes a weighting function in a weighted least square 
procedure (WLS). For every K e N we compute: 

25 

The error will be: 



(10) 

30 
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Weighting function; 

Whore A is 5 X 5 matrix. Hie rank of A is only 4. In order to 
reduce the amount of real time computation, the following computations are 
performed and the results stored 

15 (12) a=F t F 
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A =A V A 0 A? 
[AyA^eiglA) 

20 Since the rank of A is 4, one of the eigenvalues and thus the 

corresponding eigenvector are 0. Matrices 4X5 and 5X4 are created and 
designated by a ~ over the top of the matrix identifier and thus: 

(13) 

25 A=F T F 

F=A*A V T 



• 30 
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Using these values in the weighting function produces: 

(14) c(& = (35-^^(5-4 

5 

10 

The full algorithm thus becomes: 
Compute in off-line: 

1. Compute P. 

15 _ 

2 . For all K£N compute n= FK. 

3 . Compute the matrix 22*. 



Compute in real time: 

20 

1. Compute: 



25 



2. For all K e N, compute: 



C(if) = 



3. Find the best three (3) K that mosdy minimize C(K). This 
30 procedure is done for each hop. 

4. Compute X(K) for those K selected 
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5. Compute a histogram from all the X®, using for each DOA 
its distance from the lines (weight) in inverse as weight 

6. Filter the histogram by passing a boxcar low-pass filter of 3° 
over the histogram in order to cancel small peaks in the 

5 ndghboihood of a real peak 

7. Find the highest peak and the second to the highest peak 

8. Report the results. 

Triangulation Process 

10 The input to the location algorithm are reported measured 

directions from several base stations. Each base station reports two directions. 
Each direction is reported with its peak heights which is an indication of its 
validity. The location algorithm projects the direction of the incoming signal 
on a gnomonic plane tangent to earth, it filters out the false readings, and it 

1 5 computes the most probable location of the RMU. 

Cassini UTM Projection 
The solution of triangulation on the earth ellipsoid is a highly 
non-linear problem. In order to simplify the problem, we transfer the 

20 observed data to a plane called the gnomonic plane as shown in Figure 10. 
The geographic location of the base station on earth is transferred to an X-Y 
position on the gnomonic plane. Each measured bearing is also transferred to 
the gnomonic plane with a small adjustment. The estimation process is done 
on the plane. The estimated location of the remote mobile unit is not 

25 transferred back to an earth ellipsoid since the same plane is used for graphic 
display of maps using a Microsoft Windows program or the like. 

The Cassini UTM projection is widely used for mapping. Its 
advantages are that the mapping are in the Cartesian coordinate system such 
that the E-W and N-S lines are orthogonal. The projection maintains 

30 approximately both the directions related to the north and the distances at 
radius less than 50 M. The projection is very convenient for graphic display 
such as in the Microsoft Windows format. 
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The main distortion of the TJTM projection is caused by the 
convergence of the meridian toward the north There are two corrections that 
have to be made: 

5 1 . Correction of the map north to the geographical north, which is 

a bias correction; and 
2. Small correction of the DOA (direction of arrival of the signal) 
according to the location of the base station 

10 In actual observations, the correction was actually measured 

together with the calibration and measurement of the pointing of the direction 
finding array toward the north The calibration was performed relative to 
known reference point locations in the covered area. To inprove system 
accuracy, a reference transmitter having a known exact location may be used 

15 to correct the system parameters. 

Direction Finding Fixing on a Plane 
The fixing process input are the two DOA received from all the 
base stations. There are several tests performed on the results received from 
20 the BS: 

If the peaks are too low, both measurements are ignored 

If the ratio between peaks is larger than 4, the second peak is 

ignored 

The CS (Central System) computes the fix by using all the combinations from 
25 the DOA received from the base stations. Note that within one message there 
could be several DOAs. 

The fixing process is an LMS (least mean square) process in 
which the most probably point is found The most probable point is defined 
as the point on the plane that its sum of the squared distances from all the 
30 DOAs is minimal. The process is as follows: 

The distance of a point (X,Y) from a DOA is: 
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(15) 



d d (X, Y) = ( Y-Yj cos (Oi) - \X-Xj) sin (6J 



where (X;, Y;) are the coordinates of the base station number i. 
dj (X , Y) - distance of (X , Y) from the DOAnumber i. 
5 Q - The angje of DOA number i from the north 

(X , Y) - a point on the plane 

(15) could be written as: 

10 (16) XjsinWi) -Yi cos (0^ sin (Gj) -Y cosOJ +d ± 

Assuming that M inputs from the base station, (i e {1, ,M}), 

equation (16) is written in a matrix form; 



15 (17) 



X ± sin (8^ -Y ± cos (e 4 )' 



X M sinfe^) -Y M cos (8 tf ) 



sin (6^ -cos (8^ 



sinfB^) -cos (8 tf ) 



ft 



20 



where M - number of bases, as shown in Figure 9. 
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(18) 



10 



15 



X M sin (0^) -Y M cos (d M ) 
sin (0j) -cos (Si) 



sirite^) -cos(0 M ) 



2^ 

■4, 



Where: 

z is a vector of M elements. 
20 A is a matrix of M rows, 2 columns. 

c is a vector of M elements. 
D is a vector of M elements. 
The purpose is to find z such that D 1 !) will be minimized. The 
algorithm is LMS, as described above, and is performed by differentiating 
25 D T D relative to c and comparing the result to 0. 
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Complete Fixing Al gorithm 

1. Compute all intersections between all probable DOAs using 
(19). 

2. Perform two dimensional filtering over all intersections, taking 
into account the power of each intersection according to the 
DOAs power that created that intersection The filter size is 
900 X 900 metere. 

3. Find the an intersection from the DOA that are in the filter. 



10 Operation Summary 

In summary, and in referring once again to Figure 6, the slope 
of the lines is 2/3 and the algorithm finds the closest line to the observed 
point 600. The most probable actual point lies on the line 601. Figure 7 
shows the use of more ambiguity lines to get better resolution of the observed 

15 point 701. The actual probable location of the observed point would be on 
die line at position 702 where the slope of each of the ambiguity lines is the 
distance between the antennas. Thus a closer spacing of the antennas shown 
at Figure 7 produces less ambiguity as to the observed versus actual angle of 
the DOA 

20 The ambiguity lines shown in Figures 6 and 7 are for three 

antennas. In reality according to the matrices described above, six antennas 
are used in a five dimensional space which is unable to be drawn in a figure. 
In a five dimensional space, the observed point is located, and the nearest 
lines in the five dimensional space are located for the most probable actual 

25 DOA for each based station By spreading out the spacing of the antennas, 
more ambiguity lines as shown in Figure 7 are used such that antenna spacing 
is best spread out to minimize ambiguity. 

In the based station, the observed point 701 representing DOA 
0 is compared to the closest three lines in a five dimensional space. A 

30 histogram is then generated based on the three lines closest to the observed 
point Then a gaussian distribution of the observed 0s is computed and only 
the actual points 702 and five dimensional space are picked which he within a 
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single standard deviation. A histogram typically produces two peaks for most 
probable observed values. Histograms that produce more than two peaks are 
filtered so that only the two strongest peaks are taken as potential values. 
These two most probable actual values of the DOA 0 are then sent to the 
5 central she for triangulation to determine the position of the transmitter. 

Those skilled in the art will readily recognize that the benefit of 
the histogram analysis may be achieved through a less confutation intensive 
operation known as K/m where the DOA 6 most often CKxnrring in the 
analysis is the most probable DOA 6. A complete histogram need not be 
10 computed if instead a statistical threshold is applied in which the number of 
occurrences of 9 exceeds a threshold value. As soon as the .threshold value 
is readied, further analysis is halted since the most probable 0 has been 
found 

At the central site, the two most probable 0 for the two 

15 histogram peaks around a plus or minus three degree range are recorded for 
two, three or more base stations. Thus T most probable points are received 
by the base station whereas n is equal to the number of stations. Thus if four 
stations are reporting, two most probable DOA 0, 16 observed points are 
received at the central station. A sparial filter is used with a 900 by 900 

20 meter space to determine which of the points cluster within the most probable 
location. Within this 900 square meter space, the most probable location is 
selected Additional computation can be made by overlaying the 900 by 900 
square meter space onto a map of an urban area and, for vehicle location, 
detemiining that the vehicle could only be on a street and not within a 

25 building. In actual observations, the actual location of the remote mobile unit 
is within a 77 square meter area. 

Thus the multipath problem is solved through the use of a 
multidimensional ambiguity space to eliminate observed DOA values which 
are least probable. The use of a plurality of arrays will produce a number of 

30 observed values within the n- dimensional ambiguity space, some of which are 
due to the actual direction of arrival and some of which are due to multipath 
reflections. Since the observed values in the preferred embodiment of the 
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present invention are the result of multiple frequencies in a frequency hopped 
spread spectrum system, observed values and most probable values will be 
calculated at various frequencies. Different frequencies observed over a 
period of time may be more or less susceptible to the iffidttpath problem. In 
5 addition, for a moving remote mobile unit, the dbppler shift will be observed 
but minimized over time since only the most probable values due to the 
strongest signals will survive a histogram analysis: Thus die histogram 
analysis should be performed ova* a number of frequency hops over a longer 
period of time until the multipath signals statistically drop out Ibus only the 
10 true direction signals will survive the histogram analysis and will be passed to 
the central station for final analysis. 

CONCLUSION 
Although specific embodiments have been illustrated and 
15 described herein, it will be appreciated by those of ordinary skill in the art 
that any arrangement which is calculated to achieve the same purpose may be 
substituted for the specific embodiment shown This application is intended 
to cover any adaptations or variations of die present invention Therefore, it 
is manifestly intended that this invention be limited only by the claims and 
20 fiie equivalents thereof 
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WHAT IS CLAIMED IS: 

1 . A mobile transmitter location system, camprisii^ 

a mobile transmitter having a unique identifier, and having a 
5 ; transmitter for sending a frequency-hopped spread spectrum carrier signal 
containing the unique identifier, 

a first base station having an array of irregularly spaced dipole 
antennas operable for receiving the frequency-hopped spread spectrum carrier 
signal and capable of determining a first direction of origin of the frequency- 
10 hopped spread spectrum carrier signal; 

a second base station also having an array of irregularly spaced 
dipole antennas operable for receiving the frequency-hopped spread spectrum 
carrier signal and capable of detennining a second direction of origin of the 
frequency-hopped spread spectrum earner signal; and 
15 a central station connected to the first base station and the 

second base station, and operable for: 

1.) receiving the first direction of origin of the 
frequency-hopped spread spectrum carrier signal from the first 
base station; 

20 2.) receiving the second direction of origin of the 

frequency-hopped spread spectrum carrier signal from the 
second base station; 

3.) triangulating the first direction with the second 
direction to determine the location 

25 

2. The mobile location system according to claim 1 wherein a reference 
station operates periodically to transmit a reference frequency-hopped spread 
spectrum carrier signal of known origin and location to calibrate the first and 
second base stations and the central station. 

30 

3. The mobile location system according to claim 1 wherein the first base 
station and the second base station both include a reference antenna such that 
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the phase difference of the frequency-hopped spread spectrum carrier signal 
between each dipole of the array of irregularly spaced dipole antennas and the 
reference antenna is determined 

5 4. The mobile location system according to claim 1 wherein the first base 
station and the second base station both operate to reduce the ambiguity of the 
direction of arrival of the incoming signal by calculating a multidimensional 
ambiguity space having lines incfiratirig the only possible angles of the 
direction of arrival of the signal and detennining the probable points of these 
10 lines which are closest to the observed value. 

5. The mobile location system according to claim 4 wherein the first base 
istation and the second base station both generate a histogram of the probable 
values and sending the two most probable values to the central station. 

15 

6. The mobile location system according to claim 1 wherein the central 
site received the most probable values from the first and second base stations 
and applies a spatial filter to only use the most probable values which lie 
within a preselected distance of the average location of all most probable 

20 values. 



7. A method of locating a mobile transmitter location, comprising 

the steps of (", 
receiving at a first base station a radio frequency signal on a 
25 plurality of irregularly spaced antennas; 

receiving at a second base station a radio frequency signal on a 
plurality of irregularly spaced antennas; 

calculating the observed angle of arrival of the radio frequency 
signal by comparing the phase difference of the radio frequency signal 
30 between pairs of the plurality of irregularly spaced antennas; 
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determining the most probable angle of arrival of the radio 
frequency signal by comparing the observed angle of arrival to ambiguity 
lines in an n-dimensional space where the number of antennas is n+1; 

calculating an histogram of the most probable values and 
5 selecting the strongest most probable values; 

applying a spatial filter to the strongest most probable values 
from each base station and selecting the location based on the closest 
grouping of the strongest most probable values. 
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